Decoherence of orbital angular momentum entanglement in a turbulent atmosphere 
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The evolution of an entangled photon state propagating through a turbulent atmosphere is for- 
mulated in terms of a set of coupled first order differential equations, by using an infinitesimal 
propagation approach. The orbital angular momentum (OAM) basis is used to described the den- 
sity matrix of the state. Although the analysis is done in the paraxial limit for a monochromatic 
optical field, the formalism is comprehensive in the sense that it does not require any assumptions 
about the strength of the turbulence and it can incorporate any spectral model for the turbulence. 
As a comparative example the case of entangled qubit OAM biphoton states is considered. 
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Orbital angular momentum (OAM) as a basis for en- 
tangled photon states recently became popular due to the 
potential for the higher dimensional quantum informa- 
tion processing and communication that it enables [H-[3] ■ 
Since single mode optical fibre (the workhorse in the op- 
tical fibre communication infrastructure) only supports 
modes with zero OAM, free-space is the natural choice 
for OAM based quantum communication. Unfortunately 
OAM entanglement suffers decoherence in a turbulent at- 
mosphere — the random phase modulations, induced by 
the spatial fluctuation in the refractive index, scatter the 
initial OAM states into other OAM states. 

Some studies have been done on the effect of scintilla- 
tion on classical optical beams with specific OAM modes 
There has also been some work done on quan- 
tum OAM states propagating through turbulence @ [l(| , 
however, these studies modeled the turbulence by a sin- 
gle phase-only transmission function that quantifies the 
turbulence by a single parameter, the Fried parameter r . 
Reliable calculations of the effects of turbulence on the 
propagation of light require more accurate models, which 
are usually represented by the power spectral densities of 
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the refractive index fluctuation 

The approach that is followed in this paper, is to con- 
sider the incremental change in the density operator dur- 
ing an infinitesimal propagation through the turbulent 
atmosphere and to express the result as a differential 
equation. The fluctuation in the refractive index is ex- 
tremely small compared to the average refractive index. 
As a result the distance over which the effect of the 
medium becomes visible is much longer than the cor- 
relation distance of the medium. One can therefore ig- 
nore correlations in the medium along the propagation 
distance, resulting in a Markov approximation. This 
approach is reminiscent of a Master equation approach, 
but instead of a time derivative, we obtain a derivative 
with respect to propagation distance. The analysis is re- 
stricted to monochromatic optical fields under the parax- 
ial approximation. 

When propagating through a turbulent atmosphere, 



the density operator for any photon state expressed in 
terms of an OAM basis, for example for a single photon 
given by 



P( z ) = l TO ) Pmn(z) 



(1) 



would change due to the distortion of the OAM modes. 
Here p m n(z) is the density matrix elements and |m), etc. 
denotes the OAM basis elements. (The two indices of 
the OAM modes — radial index p and azimuthal in- 
dex I — are for the sake of notational simplicity com- 
bined into one index.) Since the turbulence model is 
represented by a spectral function, it is necessary to 
express the OAM basis elements as momentum-space 
wave functions G m (K,z) = (K|m,z), where K is the 
two-dimensional spatial frequency vector on the trans- 
verse momentum-space. Formally the OAM wave func- 
tions G m (K,z) are given by two-dimensional Fourier 
transforms of the Laguerre-Gaussian modal functions. 
The explicit z-dependence is a reminder that the OAM 
modes are functions of the propagation distance. For any 
given z, the OAM modes \m, z) form a complete two- 
dimensional orthogonal basis, and, by implication, so do 
their momentum-space wave functions G m (K, z). 

Another consequence of the extremely small index fluc- 
tuation is that the wave equation for monochromatic light 
propagating through turbulence separates into two terms 
in the paraxial limit: one for free-space propagation with- 
out turbulence and another for the phase modulation in- 
duced by the index fluctuation. For propagation over 
an infinitesimal distance dz, the momentum-space wave 
function is transformed as follows 



G m (K,z) A G m (K,z) + ^[|K| 2 G m (K,z) 
-2k 2 N(K,z)*G m (K,z)] , 



(2) 



where * represents convolution, and k is the wave number 
(k ^ |k|). The two-dimensional Fourier transform of the 
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refractive index fluctuation is denned as 



N(K,z) = 



*o(k) 



1/2 



cxp(— ik z 



where k is the full three-dimensional spatial frequency 
vector, $o(k) is the three-dimensional power spectral 
density of the refractive index fluctuation, and x(k) is 
a three-dimensional random complex spectral function, 
with Afc being its coherence width in the frequency do- 
main. Since the index fluctuation is a real-valued func- 
tion, x*(k) = x(— k). Furthermore, 



(X(k 1 )x*(k 2 )> = (27rA fc ) 3 5 3 (kr - k 2 ), 



(4) 



where ^(k) denotes a Dirac-delta function in three di- 
mensions. 

It is important to note that, while G m (K, z) represents 
an orthogonal basis, the distortions caused by the scintil- 
lation process denoted in Eq. @, produce a momentum- 
space wave function that is not an element of the orthog- 
onal basis anymore. The point of this approach is to re- 
express the resulting distorted wave function in terms of 



the orthogonal OAM basis and to incorporate the coeffi- 
cients of this expansion into the density matrix elements. 

First we consider the case for a single photon. The 
resulting expressions can then be used to obtain equiva- 
lent expressions for entangled multi-photon states. The 
density matrix for a single photon, in terms of the 
momentum-space wave functions, is given by 



/OO 
|Ki) G m (Ki)p mn ( 2 ) 

d 2 ^! d 2 K 2 



xG:(K 2 ) (K 2 | 



4tt 2 4tt 2 ' 



(5) 



where we now dropped the z-dependence from the 
momentum-space wave functions. 

Applying the propagation process of Eq. @ to the 
density operator in Eq. ([5]) repeatedly and averaging over 
the index fluctuation using Eqs. (|3]) and one obtains 
an expression for the change in the density operator, 



/oo />oo 
G^(K0G n (K 2 ) k 2 <&i(K')Gp(Ki — K')G*(K 2 — K') 
-oo J — OO 



d 2 K' 



ninpq 
„2 



k G p (Kx)G* (K 2 ) / $i(K 



^ d 2 K' 



4tt 2 



-flKxf-lKaHGptK^tKa) 



4tt 2 

d 2 ^Ti d 2 K 2 
4tt 2 4tt 2 : 



(6) 



where $i(K) = &o(K,k z = 0), which follows from the 
Markov approximation. 

The first term in the square brackets in Eq. ([6]) requires 
the evaluation of integrals of the form 

W mn {K,z) = J G m (K 1 )G;(K 1 - K) (7) 

The If mn 's are then used in the remaining integrals of 
the first term to give 

f°° d 2 K 
L mnpq = fc 2 ^jMK)^* p (K)l¥„ 9 (K) (8) 

The integrals over Ki and K 2 for the second term in 
the square brackets in Eq. © represent pure orthogonal- 
ity conditions and give rise to Kronecker delta functions. 
So the second term becomes —S mtP S n ^ q LT, where 

L T = k 2 J ^{K)^. (9) 

The third term in the square brackets in Eq. ©, which 
is associated with the free-space propagation of OAM 



modes, is given by 

V mnpq = — J (|K 1 | 2 -|K 2 | 2 )G P (K 1 )G^(K 1 ) 

xG„(K 2 )G*(K 2 ) ffiffi - (10) 

It is an orthogonality condition with respect to the az- 
imuthal indices of the OAM modes, but not with respect 
to the radial indices. The latter is a result of the fact that 
the OAM modes depend on the propagation distance. 

Considering the case of an entangled biphoton, one 
still finds integrals of the forms given in Eqs. (fTHTUl). 
The resulting expression for the density matrix elements 
in the case where one of the two photons propagates 
through turbulence, while the other propagates through 
free-space without turbulence, is given by 

&zPmnpq — VmnrsPrspq VpqrsPjnnrs 

~\~L rnnrs p rS pq Lj* Pmnpq > (H) 

where repeated indices are summed over. The equation 
in Eq. has the form of a Lindblad equation where 
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the first two terms on the right-hand side represent the 
Hamiltonian term and the last two terms represent the 
dissipative terms. 

In general Eq. (fTTj) represents an infinite set of cou- 
pled first order differential equations. Even if the initial 
state contains only a few lower order modes, the turbu- 
lence will cause these modes to be coupled into all other 
modes. Subsequently the other modes will couple back 
into the original modes. Truncating the set of equations, 
one inevitably excludes part of the coupling among all the 
different modes. However, this coupling should become 
progressively smaller for higher order modes. Hence, one 
may be able to truncate the set at some point while re- 
taining the dominant inter-modal coupling. 

To compare this result with previous work Q , we con- 
sider a severely truncated case, only retaining modes of 
the lowest radial index (p = 0) and with azimuthal in- 
dices of the same magnitude I — ±q. We'll consider three 
cases where q — 1,2,3, respectively. The truncation im- 
plies that the trace of the density matrix is not equal 
to 1 anymore. One can perform a normalization on the 
truncated density matrix to ensure that its trace remains 
1. However, the reduced trace gives an indication of the 
loss of information to the higher order modes. 

The formalism allows one to include any spectral model 
<f>o(k) for the turbulence, however, for the sake of com- 
parison we neglect the effect of the inner and outer scales. 
Here we use the von Karman spectrum 
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$o(k) = 



0.033C 2 



(|k| 2 



\ll/6 : 



(12) 



where C 2 is the structure constant for the turbulence 
and kq is inversely proportional to the outer scale of the 
turbulence. The outer scale will help us to regularize the 
integrals, but will disappear from the final expressions. 
Substituting Eq. (fT2|) into Eq. ([9]), one obtains, 



0.1244C 2 k~ 
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0(1) 



(13) 



For the case under consideration V m npq(z) — 0. After 
evaluating the integrals for L mnrs one finds that, in the 
limit of large outer scale, the only nonzero terms are 



L-q,q,-q y q(z) 



(z) — L q .- q . q ^ q (z) 

= L -1-q-q-q{ Z )= L T-A q h{z) (14) 



and 



-q-q{z) = L- q - qiqiq (z) = B q h(z), 



(15) 



where A q and B q are positive constants that only depend 
on q (see Table IJ), and h(z) is the same function for all 
the terms. It contains all the dimension parameters 



h(z) 



(- 
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-3 
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(16) 



TABLE I: The numerical values for the constants A q and B q 
that appear in L 



Lrnnrs TOT Q 


= 1,2,3. 






9=1 


q = 2 


g = 3 




1.570 


2.206 


2.807 


B q 


0.03030 


0.004787 


0.001754 



where zr is the Rayleigh range (ttujq/X), ujq is the radius 
of the beam waist and A is the wavelength. 

From Eq. (|13p one can see that Lt diverge in the limit 
of large outer scale. However, the nonzero elements for 
the Lx-term in Eq. (fTTj) are the same as those in Eq. (|14l) . 
As a result all the L-r-terms cancel exactly and the outer 
scale drops out of the final expression. 

Provided that the turbulence is not too weak and that 
the beam waist is not too small, the concurrence decays 
to zero over a distance much shorter than the Rayleigh 
range, which allows one to assume that 1 + z 2 jz 2 R « 1. 
Under these circumstances one can express the integral 
of h(z) in terms of the Fried parameter tq, 



[ h(z') dz' = 0.592 ( — 
Jo \ r o 
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(17) 



where r = 0.185(A 2 /C^/z) 3/ ' 5 . Thus all the dimension 
parameters are combined into wo/Vo- 

Assuming that the initial state of the density matrix is 
the singlet Bell-state in the OAM basis, one obtains the 
following solution of the density matrix 



_ T 

Pmnpq 



l-R 2 





1 + R 2 -2R 
-2R l + R 2 

l-R 2 



, (18) 



where rap (nq) denote the row (column) indices, and 
where 



T = exp 
R = exp 



~{A q — B q ) f h(z') dz' 

Jo 

-B q [ h(z') dz' 

Jo 



(19) 
(20) 



The eigenvalues of the density matrix are T(l + i?) 2 /4, 
T(l - i?) 2 /4, T(l - i? 2 )/4 and T(l - i? 2 )/4, which are 
all positive. The trace of the density matrix is given by 
T, which is a decaying function, since A q > B q , as shown 
in Table HI We plot the trace T as a function of uiq /tq in 
Fig. Ufa) for q= 1,2,3. 

Using the normalized density matrix [by setting T = 1 
in Eq. (|18[) ]. we compute the concurrence of formation 
HE! and obtain C = (2R + R 2 - l)/2, which is plotted 
in Fig. [TJb) for q = 1,2,3 as a function of wo/ r o- 

From the curves for the trace in Fig. HJa) one can see 
that modes with higher OAM are scattered more rapidly 
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into other modes than those with lower OAM. On the 
other hand, from Fig.[ljb) we see that modes with higher 
OAM retain their entanglement for longer distances than 
those with lower OAM. These conclusions agree qualita- 
tively with previous work [e}, however, while the scatter- 
ing into other modes occurs at a scale where tq k l)q, 
similar to what was found before js}, the entanglement 
lasts for at least an order of magnitude longer, which is 
quantitatively different from what was found before 
Here the slowness of the decay in the concurrence is a 
result of the smallness of the values of the BqS, given in 
Table Q] From these results it appears that the effect of 
scattering and the implied loss of photons in the desired 
OAM modes may turn out to be a more significant chal- 
lenge for free-space quantum communication than the de- 
coherence of OAM entanglement. 

(a) 
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FIG. 1: Plots of (a) the trace of the density matrix and (b) 
the concurrence for a biphoton, initially in the singlet Bell- 
state, in terms of two OAM states with I — ±q, for q = 1, 2, 3, 
as a function of uia/ro. 

In conclusion, we derived an infinitesimal propagator 
equation, reminiscent of a Master equation, that models 
the spatial evolution of a density operator for a OAM 
entangled multiphoton state, propagating through a tur- 
bulent atmosphere. The approach assumes monochro- 
matic light in the paraxial limit, but does not require 
any simplification to the turbulence model. As a result 
this formulation is capable of analyzing the effect of re- 
alistic atmospheric turbulence on the propagation of any 



quantum photon state in the OAM basis. The resulting 
infinite set of first order differential equation was trun- 
cated to consider an example for comparison with previ- 
ous work. Although the result was found to be qualita- 
tively similar, significant quantitative differences exist. 
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